A unique spinodal region in asymmetric nuclear matter 
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Asymmetric nuclear matter at sub-saturation densities is shown to present only one type of 
instabilities. The associated order parameter is dominated by the isoscalar density and so the 
transition is of liquid-gas type. The instability goes in the direction of a restoration of the isospin 
symmetry leading to a fractionation phenomenon. These conclusions are model independent since 
they can be related to the general form of the asymmetry energy. They are illustrated using density 
functional approaches. 
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Phase transitions are universal phenomena of matter in 
interaction. The coexistence regions which correspond to 
thermodynamically-forbidden areas exhibits general fea- 
tures such as metastabilities or instabilities. Since strong 
interaction between nucleons is of Van der Waals type 
nuclear systems are expected to present a liquid and a 
gas phase characterized by their respective densities 
Since nucleons can be either protons or neutrons the tran- 
sition occurs in a two fluid systems which may lead to a 
richer phase diagram involving up to two order parame- 
ters. It has been argued that, asymmetric nuclear mat- 
ter (ANM) presents different types of instabilities : a 
broad chemical instability region with the concentration 
as order parameter and a narrower domain of mechanical 
instability for which the total density plays the role of a 
second order parameter In this letter, we will show 

that spinodal instability and phase transition in ANM 
involves a unique order parameter which reflects density 
fluctuation between the two phases. We will stress that 
the transition induces variations of the concentrations 
leading to the isospin fractionation 0| similar to the one 
experimentally observed These properties and in fact 
all the characteristics of the instabilities are independent 
of their usual classification as chemical or mechanical in- 
stabilities. These conclusions are related to general prop- 
erties such as the isospin dependence of the asymmetry 
energy. The shape of the spinodal region and the asso- 
ciated instability times depend upon the model so the 
observation of spinodal decomposition may provide con- 
straints on the isospin dependence of the effective forces. 

Let us consider ANM characterized by a proton and 
a neutron densities pi — pp, pn- These densities can be 
transformed in a set of 2 mutually commuting charges 
Pi = pi, P3 where pi is the density of baryons, pi — 
Pp + Pn, and p3 the asymmetry density pa = p„ — pp. 
In infinite matter, the extensivity of the free energy im- 
plies that it can be reduced to a free energy density : 
F{T,V,N^) = VT{T,pi). The system is stable against 
separation into two phases if the free energy of a sin- 
gle phase is lower than the free energy in all two-phases 
configurations. This stability criterium impHes that the 
free energy density is a convex function of the densities 



Pi. A local necessary condition is the positivity of the 
curvature matrix : 
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where we have introduced the chemical potentials pj = 
^^\T,v,Ni = §^\T,pijij- In the considered two-fiuids sys- 
tem, the [J^ij] is a 2 * 2 symmetric matrix, so it has 2 real 
eigenvalues j3| : 



= ^ ( Tr [T,,] ± VTr [:F.,,r - 4Dct [T,, 



associated to eigenvectors 6p defined by {i ^ j) 
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Eigenvectors associated with negative eigenvalue indicate 
the direction of the instability. It defines a local order 
parameter since it is the direction along which the phase 
separation occurs. The eigen values A define sound veloc- 
ities, c, by (? = j^Pi A. In the spinodal area, the eigen 
value A is negative, so the sound velocity, c, is purely 
imaginary and the instability time r is given by r = d/\c\ 
where d is a typical size of the density fiuctuation. 

The requirement that the local curvature is positive 
is equivalent to the requirement that both the trace 
(Tr[jrjj] = A"*" -I- A~) and the determinant {Det[Tij] = 
A+A^) are positive 



TT[Tij] > 0, and Det[J^y] > 



(4) 



The use of the trace and the determinant which are two 
basis- independent characteristics of the curvature matrix 
clearly stresses the fact that the stability analysis should 
be independent of the arbitrary choice of the thermody- 
namical quantities used to label the state e.g. {pp, pn) 
or (pi. Pa). If Eq. violated the system is in the un- 
stable region of a phase transition. Two cases are then 
possible : i) only one eigenvalue is negative and one or- 
der parameter is sufficient to describe the transition or 
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ii) both eigenvalues are negative and two independent or- 
der parameters should be considered meaning that more 
than two phases can coexist. 

For ANM below saturation density, the case ii) never 
occurs since the asymmetry energy has always positive 
curvature (^33). Indeed, the asymmetry term in the 
mass formula behave like {N — Z)"^ times a positive func- 
tion of A showing that the dominant ps dependence of 
the asymmetry potential energy is essentially quadratic 
and that .^^33 is a positive function of the total den- 
sity. Recent Bruckner calculations in ANM have con- 
firmed the positivity of .F33. They have parameterized 
the potential energy with the simple form V(pi,p3) — 

Vo{pi)pl + Vi(pi)pi with Vo/Vi 3 and Vq < 0. 

This is also true for effective forces such as Skyrme 
forces. For example, the simplest interaction with a con- 
stant attraction, io, and a repulsive part, t^pi, leads to 
V{pi,p3) = {3pl-pl)(3{pi) with p{pi) = (io + t3/6pi)/8. 
The function /3{pi) is negative below saturation density, 
hence the contribution of the interaction to in the 
low density region is always positive. 

These arguments show that, below saturation density, 
the p3 curvature, ^^33, is expected to be positive for all 
asymmetries. Since the curvature in any direction, Tu, 
should be between the two eigenvalues A" < < A+ we 
immediately see that if JF33 is positive one eigen curvature 
at least should remain positive. In fact for all models 
we have studied appears to be always large enough 
so that the trace is always positive demonstrating that 
A"*" > 0. Since Tr[J^ij] = JF„„ -|- J-pp, this can be related 
to the positivity of the Landau parameter jr„„ and Tpp. 

The large positive value of T33 also indicates that the 
instability should remain far from the direction i.e. it 
should involve total density variation and indeed we will 
see that in all models and for all asymmetry the insta- 
bility direction hardly deviates from a constant asymme- 
try direction {Sp^ ^ Spi). This isoscalar nature of the 
instability can be understood by looking at the expres- 
sion of the eigen- modes in the (p„; pp) coordinates. Since 
A~ < J^ii < A"*", the differences A" — J-u is always nega- 
tive demonstrating, using Eq.|Sl that the instability is of 
isoscalar type, if the Landau parameter Tnp is negative. 
Hence, there is a close Hnk between the isoscalar nature 
of the instability and the attraction of the proton-neutron 
interaction 0. 

To illustrate the above results we will now use density 
functional formalism using Skyrme and Gogny [lfli | ef- 
fective forces. It should be noticed, that the extraction 
of the sound velocity corresponds to a random phase ap- 
proximation (RPA) which goes beyond the mean field 
approximation. We will focus on the isospin degree of 
freedom for which RPA approaches can be considered 
as good approximations and discuss both zero and finite 
temperature in ANM. 

We represent in Fig.^the energy surface as a function 
of Pn and Pp, deduced from SLy230a Skyrme interac- 
tion In the symmetric case (p„ = Pp), one can see 
the negative curvature of the energy which defines the 
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FIG. 1: This figure represents the energy surface as a function 
of the densities p„ and Pp for the SLy230a interaction. The 
contour delimitate the spinodal area. 



spinodal area, whereas in pure neutron matter {pp = 0), 
no negative curvature and so no spinodal instability are 
predicted. We can also notice that the isovector density 
dependence is almost parabolic illustrating the positivity 

of .F33. 




FIG. 2: This two figures are a projection of the spinodal con- 
tour in the density plane : left, for Skyrme (SLy230a [lljl. 
SGII J0, SIII [11) and Gogny models (Dl 0|, DIS IH, 
DIP [1^) ; right, temperature dependence of the spinodal 
zone computed for the SLy230a case. 

The spinodal contours predicted by several models ex- 
hibit important differences (see Fig. EJ. In the case of 
SLy230a force (as well as SGII, DIP), the total den- 
sity at which spinodal instability appears decreases when 
the asymmetry increases whereas for SIII (as well as Dl, 
DIS) it increases up to large asymmetry and finally de- 
creases. Considering the impHcit equation of the spinodal 
limit g{pi,P3) = Det[!Fij] = 0, we can show that the cur- 
vature of the spinodal around the symmetry is related 
to d^g/dp^{p3 = 0). The isospin symmetry impose that 
dg/dpz{pz = 0) = and one can show that 



5^1 
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IP3=0 



(•^1133-^33 + •^3333-^11 + •?^1333-^13)- (5) 
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Assuming that the forth derivatives of the free energy 
are negUgible compare to the third one, the curvature 
of g is mainly given by T>e%[dTij / dp^] where dTijjdpz 
is nothing but the curvature matrix of dF jdpz., i.e. the 
asymmetry density dependence of the free energy. We ob- 
serve that all forces which fulfill the global requirement 
that they reproduce symmetric nuclear matter (SNM) 
equation of state as well as the pure neutron matter cal- 
culations, leads to the same curvature of the spinodal 
region. 

The temperature dependence of the spinodal contour 
can be appreciated in the right part of Fig. ISJ. As the 
temperature increases the spinodal region shrinks up to 
the critical temperature for which it is reduced to SNM 
critical point. However, up to a rather high temperature 
( 5 to 10 MeV) the spinodal zone remains almost identical 
to the zero temperature one. 
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FIG. 4: This is the projection of the iso-eigen values on the 
density plane for Slya (left) and DIP (right). The arrows in- 
dicate the direction of instability. The mechanical instability 
is also indicated (dotted line). 
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FIG. 3: Sound velocity as a function of the density. Negative 
values of sound velocity means that it is imaginary (e = i is 
c~ imaginary). On the left, we have changed the asymmetry 
parameter and fixed the model (SLy230a), on the right, we 
have changed the models and fixed y = pp/pi=0.5. 

We show in the right part of Fig.Elthe sound velocity 
in SNM as a function of the density for several forces. 
When we enter into the spinodal area, the sound velocity 
becomes purely imaginary. The different forces predict 
different instability time. However, for the set of forces 
fitted to reproduce the symmetric and the pure neutron 
matter calculation, we observe a convergence of predic- 
tion as we already observed on the spinodal boundaries. 
In the left part of Fig. we can appreciate the reduction 
of the instability when we go away from SNM. However, 
large asymmetries are needed to induce a sizable effect. 

Various contour of equal imaginary sound velocity are 
represented in Fig.|4lfor SLy230b and DIP interactions. 
The more internal curve correspond to the sound velocity 
i0.09c, after comes ?0.06c, iO.OSc and finally 0, the spin- 
odal border. For these two recent forces which takes into 
account pure neutron matter constraints, the predicted 
instability domains are rather similar. We observe that 
in almost all the spinodal region the sound velocity is 
larger than 0.06c. 



Let us now focus on the direction of the instability. If 
bp~ is along y = pp//9i=const then the instability does 
not change the proton fraction. For symmetry reasons 
pure isoscalar (5p3 — 0) and isovector (Sp\ — 0) modes 
appears only for SNM so it is interesting to introduce 
a generalization of isoscalar-like and isovector-like modes 
by considering if the protons and neutrons move in phase 
(6p~5pp > 0) or out of phase {Sp~Sp~ < 0). Fig.01shows 
the direction of instabilities along the spinodal border 
and some iso-instability lines. We observed that insta- 
bility is always almost along the pi axis meaning that it 
is dominated by total density fiuctuations even for large 
asymmetries. Fig. presents the angle of the eigen state 
Sp~ with the isoscalar axis normaHzed by the angle be- 
tween the y=const line and the isoscalar axis (denoted x) 
for 3 models (DIP, SGII, SLy230a). We can see that this 
quantity is in between and 1, so that the instability di- 
rection is between the ?/=const fine and the pi direction. 
This shows that the unstable direction is of isoscalar na- 
ture as expected from the attractive interaction between 
proton-neutron. The total density is therefore the dom- 
inant contribution to the order parameter showing that 
the transition is between two phases having different den- 
sities (i.e. Hquid-gas phase transition). The angle with 
the pi axis is almost constant along a constant y line. 
This means that as the matter enters in the spinodal 
zone and then dives into it, there are no dramatic change 
in the instability direction which remains essentially a 
density fiuctuation. Moreover, the unstable eigenvector 
drives the dense phase (i.e. the liquid) towards a more 
symmetric point in the density plane. By particle con- 
servation, the gas phase will be more asymmetric leading 
to the fractionation phenomenon. Those results are in 
agreement with recent calculation for ANM Q and nu- 
clei 113. 

We want to stress that those qualitative conclusions are 
very robust and have been reached for all the Skyrme and 
Gogny forces we have tested (SGII, SkM*, RATE, Dl, 
DIS, DIP,...) including the most recent one (SLy230a, 
DIP) as well as the original one (like SIII, Dl). 
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FIG. 5: Angle of Sp~ with the isoscalar axis normalized by 
the angle between the j/=const line and the isoscalar axis 
(denoted x) for various values of y noted on the figure and 
for 3 models (solid: SLy230a, dotted: DIP, dashed: SGII). 

A different discussion can be found in the litera- 
ture Therein, a chemical and mechanical stability 
conditions 

^|t,p>0 and I^It., >0 (6) 
oy opi 

are introduced and two types of spinodal regions are de- 
fined. This leads to the idea that two order parameters 
should be introduced : the concentration in the chemical 
instability zone and the baryon density in the mechan- 
ical instability region. To try to connect the eigenstate 
analysis with this discussion one can use the relations 

p„Det[^,,] = ^hp|£k, (7) 

Comparing expression 0, with the identity Det[J^y] = 
A+A^ one can be tempted to relate separately ^^|t,p 
and ■§^\T,y to the two eigenvalues A"*" and A~. This is 
indeed correct in SNM Q but theses relations break down 
in ANM. For instance, ^f^lr.j/ = Pi^^It,^ is nothing but 
the curvature of the free energy in the particular direction 
of constant proton fraction and this direction has no rea- 
son to be an eigenvector direction, except in SNM. Con- 
sidering Eq.Qonly the product ^^|t,p * 'Wpll'T^y should 



be used to spot the instability region. Eq.d shows that 
at the onset of instability where r)et[Tij] = the product 
should vanish. On the spinodal border, §^\T,y vanishes 
only if the direction of instability is the ?/=const line. 
This happens only for SNM therefore in general ^^|t,p 
should vanish first irrespectively of the actual nature of 
the instability. In particular, the previous eigen vector 
analysis shows that in the so called chemical region not 
only the concentration fluctuations, Sy, are amplifled but 
mainly the total density one, 5pi. When entering in the 
mechanical spinodal zone (shown in Fig. 01 nothing spe- 
cial happens : the instability strength evolves smoothly 
(see also Fig.EJ and the associated vector keeps pointing 
in the same isoscalar-like direction (see Fig. [^. Conse- 
quently, the chemical and mechanical stability conditions 
should not be considered separately but combined into 
the determinant of the curvature matrix. 

In this paper, we have shown that ANM does not 
present two types of spinodal instabilities, a mechani- 
cal and chemical, but only one which is dominantly of 
isoscalar nature as a consequence of the negativity of the 
Landau parameter J-np. This general property can be 
linked to the positivity of the symmetry energy curva- 
ture J-33. This means that the instability is always dom- 
inated by density fluctuations and so can be interpreted 
as a liquid-gas separation. The instabilities tend to re- 
store the isospin symmetry for the dense phase (liquid) 
leading to the fractionation of ANM. We have shown that 
changing the asymmetry up to pp < 3pn does not change 
quantitatively the density at which instability appears, 
neither the imaginary sound velocity compared to those 
obtained in SNM. All the above results are not qualita- 
tively modifled by the temperature which mainly intro- 
duce a reduction of the spinodal region up to the SNM 
critical point where it vanishes. The quantitative predic- 
tions concerning the shape of the spinodal zone as well as 
the instability times depends upon the chosen interaction 
but converge for the various forces already constrained to 
reproduce the pure neutron matter calculation. 
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